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1 . Introduction 


The aim of this research is the application of the indicial re- 
sponse approach, a modeling approach originally used for studying 
nonlinear problems in flight dynamics and exposed in a recent re- 
view article [1] , to the study of vortex- induced oscillations phe- 
nomena. The problem chosen, vortex shedding from a bluff body, is 
itself a subject with many practical applications and has attracted 
the attention of researchers for over several decades. In spite 
of the long intense interest, researchers in this field agree that 
an adequate rationally based mathematical model of the phenomena 
still does not exist [2] . The phenomenon of vortex shedding en- 
ters into a class of problems involving a periodic equilibrium state 
which is the result of a bifurcation from a previous steady equi- 
librium state as a parameter exceeds some critical value. The oc- 
currence is called a Hopf bifurcation. In [1] it was shown on the 
basis of physical reasoning that the amplitude and frequency of the 
periodic equilibrium state are determined by time invariant con- 
ditions and it was suggested that specification of its phase de- 
mands a proper consideration of fluid dynamics effects. 

To that end, this study has two purposes. The first purpose is to 
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develop the indicial response approach appropriate for the prob- 
lem considered. The indicial response approach, which takes its 
origin from linear theory, enables one to write the response of a 
(linear) system as a superposition integral of step responses. This 
approach transposed to nonlinear systems in flight dynamics gives 
an analogous result in terms of a generalized superposition inte- 
gral of elementary indicial responses [1] . We give in this report 
a derivation of the above result in a suitable form for the prob- 
lem considered, e.g. in terms of the appropriate physical quan- 
tity which is the velocity field. Such a derivation provides us 
the opportunity of analyzing the mathematical postulates involved, 
which until now haven’t been properly discussed. We demonstrate 
that this main theoretical result can be applied first to the sim- 
ple situation of a time -invariant equilibrium state. The indicial 
response of the velocity field is derived directly from the equa- 
tions governing the fluid motion. These are taken to be the in- 
compressible Mavier-Stokes equations. Results for the aerodynamic 
response are shown to confirm the form, obtained previously by a 
variety of approaches [1], The analysis is then directed to the 


new situation where the equilibrium state is periodic in time. The 



second purpose is to apply the results of the analysis to the par- 
ticular case of flow past a cylinder in periodically forced motion 
to derive some salient features . The indicial response approach 
is shown to be equivalent to the currently popular amplitude equa- 
tion approach. 

2 . Generalized superposition integral based on indicial response 
approach . 

We derive the indicial response approach with emphasis on the phys- 
ical postulates involved. 

Poj- generality, consider an aircraft that has started from rest in 
the distant past with fixed axial velocity U () and zero vertical ve- 
locity ( cf . Fig. 1). Its motion is referred to an X,Y coor- 
dinate system that is fixed in space. It passes through the ori- 
gin at the arbitrarily chosen initial instant £ = 0, maintaining 
the constant axial velocity U 0 and simultaneously translating ver- 
tically, with vertical velocity v c at the center of gravity being 
an arbitrary function of time £ . The angle-of -attack a is defined 
as the angle between the resultant velocity vector and the aircraft s 
longitudinal axis: 
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( 1 ) 


a — tan 


-i MO i 
t/o 1 


Let us note that we specify a constant axial velocity U 0 to be in 
accord with normal operating conditions in wind- or water-tunnel 
experiments, wherein the uniform oncoming flow, normally held at 
constant velocity, would supply the corresponding value of Uq. To 
form the indicial response, we need to consider two motions. In 
the first one, the aircraft undergoes a variation of angle of at- 
tack a(£) from time zero to a time f = r ( cf . Fig. 2). Sub- 
sequent to time t, the angle of attack is held constant at ce(r). 
The first motion history therefore is designated as « ,(0: 


\ 




if 0 < £ < t 


( 2 ) 


[ = a(r) : if £ > t 

Notice that the subscript r is used to distinguish the above mo- 
tion history from the natural motion history a(£). In the second 
motion, the aircraft undergoes the same angle of attack history 
up to time r. Subsequent to r, the angle of attack again is held 
constant, but is given an incremental step change Ace over its pre- 
vious value of a(r) . The second motion, designated a’(£), is rep- 



resented as 


°£(0 = tt r (0 + C .»? r (0 


( 3 ) 


with : 


f = A a 


*?r(0 = 


0 : if 0 


^ < r 


( 4 ) 

( 5 ) 


[ 1 : if £>t 

Let us define by u{x, t,r)(.x : spatial coordinates, t\ time of ob- 
servation subsequent to r) the velocity field in the neighborhood 
of the aircraft. 

The indie ial response approach is established by making use of the 
following postulates. 

Post u lat e §1 : Corresponding to a motion history belonging to 
the family of motion histories « r (£) with r € |0, /]. there exists 
a velocity field u{x,t,r) well defined and supposed known by some 
means . 

In mathematical terms, one could say that there exists a mapping 
of the space of scalar quantities a T (£) with r 6 [0, t\ to the space 
of vectors u(x.t,r). However, in physical terms used in this re- 
port. we simply say that there exists a velocity response u to a 
motion history a r . An analogous postulate is held about the ex- 



istence of velocity fields associated with the family of motion his- 
tories «*(£) with r <E [0, < | : 

. Postulate §2: To each member of the family of motion histories 

a*(£) with r€ [(),(], there exists a velocity response u[x, t, r] . 

Such velocity response reproduces the velocity response correspond- 
ing to the motion history ct r+AT (£) within a negligible error of i?(Aa 2 ) : 

u[a T (f) 4 - Aa. 7 ? r (f)] = tt[o T+Ar ({)] + t?(Aa 2 ) t = l,...,n ( 6 ) 

The velocity of the angle-of -attack 6t(r) = ^ is supposed to be 
defined for each value of r 6 [0,/]. A further postulate is made 
about the existence of the indicial response of the velocity field: 

Pos tulate § 3: For every value of [0.£], there exists an in- 

dicial response u Q |x, a(f); t, r\ defined as: 

ujx, a(f);<,ri = h'm Aa - 07- {ujx, a r (f) T Ac*.77 T (£); t, t \ - u|x, a T (C); r l} 

UO! 

( 7 ) 

Given a motion history a(£ : £ 6 -ocj ), it is possible to choose 
an initial time and make a time partition * ■ ■ ,t n \ of the time 

interval [t (h t ] such that: 

t n — t , nAt = t . t r t x _ ] -= At n) 

7 


( 8 ) 



One has the following relations: 


unatj^) + Aa.i 7 t „(£) j — u Q foc t „(^)j = u a |a t „(0]Aas. + t?(Aa 2 ) 
“»i a ti(0 + Aa.i/!,(0] - «i[«t|(0] = «ai a ti(f)]Ao' + ^(Aa 2 ) 

«ni«i„- J (0] + Aa.ry (fi _, (£)j - u n _i \a tit _ l (£)] = «*[<**„_, (01 Aa + ^(Ao 2 ) 

( 9 ) 

Notice that the additional dependence of u t on x,t,t,(i = l,...,n) 
is omitted for reason of abreviation. 

By summing up the relations (9) and using postulate §2, one gets 
“*[<*<,.(£)! - Uo!<*<„(£)] = 51 “»|«e.(0]Aa + i?(Aft 2 ) (10) 


i . e . : 

u|x*a^(^); t,t I 


C ^ dot 

u (1 !x,a t „(^);/,/„ \- dr u a [x,a T (Z)-,t,T\ — + tf( Aa 2 ) (11) 
Jt„ UT 


The relation (11) constitutes our main result that will be used in 
the subsequent sections. Notice that it relies on three physical 
postulates which have to be satisfied. 


3 . Time- invariant equilibrium state 
We wish to derive the aerodynamic response on the basis of rela- 
tion (11). In the following subsections, we shall demonstrate, on 
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the basis of the Navier-Stokes equations, that it is possible to 
find analytical expressions of the velocity responses over the time 
interval t — r>0 to the motion histories «,-(£) and a*(£). 

3.1 Construction of the velocity response for the motion history (*,.(£) 
For simplicity, we neglect compressibility effects and assume that 
the fluid motion is governed by the Navier-Stokes equations for an 
incompressible fluid. In a coordinate system attached to the body, 
the Navier-Stokes equations have an additional term v c (where vc 
is the vertical velocity at the mass center) to account for the ac- 
celeration of the coordinate system relative to inertial space: 

f)lt _ _ 

— + (u.V)u + Vp-i/V 2 u = -v c (12) 

dt 

Here, u is the velocity field, p the pressure and u the kinematic 
viscosity. We wish to derive a form for the velocity response to 
the motion a r (£) over the time interval / r > 0. Since a T (£) - 
a(r) = constant for t — T> 0, boundary conditions determining the 
velocity response over the interval t — T > 0 are perfectly steady. 
Consequently, we expect the velocity response to attain an equi- 
librium state as the interval becomes large, i.e., as t—T — > oo . 

The response of the velocity field can be decomposed into an equi- 
librium state VL equt i and a transient component u, rtmJ that decays as 
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time increases. The principal condition that we impose in this sec- 
tion is that the equilibrium state u e<JU ii be time- invariant . Thus, 
as time is referred to t (the instant specifying when the angle- 
of-attack is to be kept constant thereafter), we have: 


ll(x,t T 0) Hequi ^k V'transi.Xyt T 0) 

p{x. t T > 0) — Pequil (^) ~k Ptransi.^'it T 0) 

To simplify the notation, we define: 


(13) 


t + — t — T > 0 


(14) 


Substituting (13) into (12), we get for the equations governing u (ra „.. 


""k ( u equil ■^) u tran.i "k { u trans ■^) u equil "k ^ Ptranf * u trana 

+ (u fran ,,. V)u( rans + Vf ( r ) . ( 1 — H(t + )) ~ 0 | 


^ ‘H'tranfi ® 


(15) 


where //(<_) is the Heaviside function: 


H[U) = 


0 : *+ < 0 

1 : 4 > 0 


(16) 


with utmn* = 0 on appropriate boundaries dCl. 

To find a form for u (ran? , we shall first consider a small distur - 
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bance u which satisfies the linearized version of (15) : 

37 -u + (u fi 9 Uti .V)u + {u.V)u eqml + Vp - i/V 2 u - ( r ) . ( 1 - H(t 4 )) 

V.u 1 — 0 

(17) 

with u = 0 on <9fl . Equations (17) constitute a linear eigenvalue 
problem. Since u equt i is independent of time, the eigensolutions 7 n 
are also independent of time. It is possible to choose them such 
that [3] 

u = u n - exp(X n t + )^ n {x) , V.^ n = 0 , -fnlan = 0 (18) 

From equation (18), the corresponding expression for p is: 

P = exp(A„4)p^(x) (19) 


The eigensolutions 7 n (x)(rc = 1, . . . , N) are associated with the eigen- 
values A n , the latter all having negative real parts in the phys- 
ical situation under consideration. There is an equation adjoint 
to (17) having a set of eigensolutions 7 * such that [3] 

= 6ij 7 i.J 

(liilj) = 0 V i, j 

where the brackets ( ) denote the scalar product over space. The 

eigensolutions 7 „(x) span a complete functional space and one can 

1 1 




use this fact to construct a suitable solution form for u trans (x. ) . 
Returning now to the full nonlinear equation (15) governing u tran ? . 
let us assign u tranf the form resulting from its projection onto the 

functional space of the : 

N _ _ 

Utrans{x,t + ) = Y ( d n {t + )^ n {x ) + d n {t+)^ n [x)) (2l) 

n— 1 

where the barred terms denote complex conjugates (c.c.). It will 
be convenient to break the pressure term p tra „., in (13) into two parts: 
one following the form of equation (21), and the other (denoted p) 
corresponding to the nonlinear part of the disturbance equations: 

p tran ,(x,t + ) = Y [dn{t+)p' n {S) +d n {t + )p n {x)) -T p{x,t + ) (22) 

n— 1 

After inserting (21) and (22) into (15) and eliminating the set of 
terms that identically satisfies the linear form (17), we find that 
the coefficients d n [t + ) satisfy the following equation: 

Y ^ n d n ) *)„ + — \„d n )~f n 4- Vp — — [d n d m (^i n .V)~j m + 

n n n m 

d n d m (i n .V) 7 m + d n d m tf n .V)7 m + d n d m (\.V)^ m - vr(r).(l - H(t + ))) (23) 

Multiplying this equation by the vector 7*(x). integrating over space 
and using the properties of adjoint vectors defined by (20), we ob- 
tain 

dj[t ^) Xjdj = ^ ^ ^ {^d n d m (('7n-^)0m ) ) "t d n d m ( (qf n .V) 'y Tn , 

n m 
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( 24 ) 


d'n^m \ ( ('"Yn *^) i Yji / ) 

-c,(a(T )).(]-//((.)) J- 1.....A' 
with : 

<*(««) = f f* ?j -rh, - < 25 > 

J V II v c II 

Notice that the pressure gradient term Vp in (23) makes no appear- 
ance in (24) . The properties of 7 „ and 7 * (solenoidal and vanish- 
ing on SQ) ensure that scalar products of gradient terms (such as 
Vp) and 7 * will be identically zero. 

By defining appropriate coefficients 4y nffll ^jnm^ (^jnm and Djnm* the 
system of equations (24) can be rewritten as: 

dj{t + ) - A J -(a(r))rf i (< + ) = H Minm(a(7-))d„(*-) d m(^ ) + 

n m 

Bjnm( fi { T ))d-n(t+)d m (t~) + C ; nm ( o ( t ) ) d n (t + ) d m (t j. ) -r •£^) nm (o(r))d n (£ + )d m (t + ^ 

- C> (a(r)).(l-J 5 r(< + )), j = h...,N (26) 

Thus, it appears from the system of equations (26) that it is re- 
quired to fix the values of dj(t+ = 0 )(j = 1 ,...,A T ) at some ori- 
gin of time. 

The values of dj(t + = 0) must be determined from a match of the ve- 
locity fields on either side of £ = r : 

u(x. £ = r_) - u(f, s ~ T + ) (27) 
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On the negative side of £ = t , in general the velocitj r field is 
the resultant of the entire past history of the motion a(f) up to 
the ‘‘present’’ time £ = t . Through specification of the d } at 
t + = 0 in terms of this velocity field, the velocity response for 
>0 acknowledges its dependence on the past motion. The rep- 
resentation of this past motion in a suitable way is all that re- 
mains to be done to realize a complete characterization of the ve- 
locity response for t + > 0 and thus to satisfy postulate §1 . We 
shall defer a discussion of this important step until we have de- 
rived the velocity response to the motion a*(f) and formed the in- 
dicial response. For the present, we simply designate the depen 
dence of dj(t. = 0) on the past motion as a functional: 

dj[t, =0) «(£:£€ [0,r]),* + = Oj (28) 

3.2 Construction of the velocity resp onse for the motion history <**(£) 
The velocity response fo the motion history o:*(£) can be constructed 
in an analogous way to that in the previous section. It is pos- 
sible to demonstrate that the velocity response to the motion his- 
tory a*(£) reproduces the velocity response to the motion history 
Or+Ar(0 within some negligible error of i?(Aa 2 ) for time t > t 4- 
Ar. Postulate §2 is therefore satisfied. We refer for complete 
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details to reference [4] . It is shown that the velocity response 
to the motion history or* ( ^ ) is; 


u|z,a;(£)! = u cgu n(x,a(T) - Aa) + X! d 2nK(0><+r7n( f > a ( r ) ~ Aa) 

n 


+ C.C . 


(29) 


3.3 Indicial and total responses of the velocity field 
The indicial response of the velocity field can be obtained as 


tt[£, <*;(£)] - u(f,Q r (0] 

u a = lim — — — — - 

Aa— ►(') A Ot 


(30) 


Indeed, it is permissible to make the operation defined by the right - 
hand-side member of equation (30) on the analytical expressions ob- 
tained for u[a T (f)] and uja*(f)] in the previous subsections. There- 
fore, we have satisfied postulate §3. The operation in equation 
(30) leads to the following expression for u a : 

,. v egui i(x. a(r) - Aa) - u eifUll (x,a(r)) 1 

u a — lim — — — — r lim — — . 

A Q Aa— 0 A a 

Y ( d 2r»[«2(0]'7n(f,a(’') + Aa) “ d n ; «l(0]Tn(^, <*(0)) + C.C. (31) 

n 

One gets: 


d 3 \ 

Y ( rf 2n|«2(0 - rfn|Ol(0])5n + C.C. 


(32) 
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The total response of the velocity field is obtained by summing the 
indicial response along the motion history a ( (£) according to equa- 
tion ( 11 ) : 


«(x,t) = <u.i(x,a( 0 ) + £^^j^^ n [a T (e)]^Tn(x,o : (r)) 

+ T, I'm 

^ Aa-0 


< ^2n} Q! T (01 a r(Q] 

A a 


+ c.c. 


(33) 


It has the following analytical form 

u{x,t) = UequAxi <*(«)) + J q dr^ ? \d, \a T (£) ,t + = 0 ];i + ] (34) 

with the functional dependence of J determined by the initial con- 
dition dj(t * = 0)0 = 1,. ..,N). 

3.4 Determination of the initial condition dj[t + — 0)(j — N) 

As it appears in equation (33), the value of the velocity field con- 
tains terms appearing under an integration over time such as 


exp{A > (a(r))(< - r)} 

due to the behavior of dj\a T (£)}(j = l,...,n) (cf . [4]). Hence, only 

0 j* 0 C 0 nt past motion history is taken into account in determining 
the value of the velocity field. Under such circumstances, one can 
consider a fit of the recent past motion history to a known motion 
history of the system. A motion history of the system is provided, 
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for instance, by making experiments under external periodic forc- 
ing . It can be shown, either theoretically or experimentally, that 
the response of the periodically forced system is periodic in time 
with the same period u as for the external forcing: 

u(f,r_) = ]Ta n e mT - 

n 

- a 1 (x)e ,,/T - = u(f,a(r),a(r)) (35) 

The choice of d 3 {t + =0) is governed by the following equation, de- 
rived from equation (27) 

u(x,r.) = u equtl {x,T+) + ^d„(£ = r+)'7 n (f,or(r + )) (36) 

n 

i.e. by using the properties of adjoint vectors: 

dj(£ = r + ) = (tt(f.r_),' 7 *(x,a(r + ))) - (u e<}ull (x, r+), 7 *(x, a(r + ))) (37) 

The analytical dependence of dj(t+ =0) is therefore: 

d } (t + = 0) ~ </>(a(r),a(r)) (38) 

For slow motion where o(r) ~o(7‘), one can expand d : : 

dd 

dj{t+ = 0) = d>(a(r),0) + + ,? (« 2 ( r )) ( 39 ) 

3 . 5 Determination of the lift force ~ Discussions 
Using relation (39) valid for slow motion, one can rewrite rela- 
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tion (34) as: 


u(x,t) ~ u e<Jtl ,i(f,a(f)) + f dT~I{a{r),a(T),t + ) (40) 

Jo ar 

It is possible to derive the pressure from the velocity field by 
using the Havier-Stokes equations. The value of the lift force L 
is obtained by taking the value of the pressure on the body sur- 
face ( dB ) and by integrating over a solid angle: 

L = [ p(x,t)\ 3B s\n(0)d6 (41) 

J — 7T 

Hence, the lift coefficient is of the following form 

C L (t) =C’r'\a(t )\- = 0| ,*+,rj (42) 

where the value of the function 7 approaches zero as t — t o o, 

due to the properties of d n (t^) . The relation (42) was derived in 
previous work on the basis of functional expansion approach [1] : 
the function 7 was called the deficiency function. It was used 
for various applications in aircraft dynamics (cf. for instance 
[5]). The value of the lift force can be approximated as: 

C L (i) ~ - J o dT^?(a(T),a{r),t+) 

- Cr i M01-/ o ‘*^/(a(r),4) + ^(a 2 (0) (43) 
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Such relation is used in most mathematical models in flight dynam- 


ics . 

4. Periodic equilibrium state 

As one of the values of the real part of the eigenvalues A ; (j = 
l,...,A r ) becomes positive for increasing Reynolds number, the time- 
invariant equilibrium state loses its stability. A Hopf bifurca- 
tion occurs: the time-invariant equilibrium state is replaced by 

a time-varying equilibrium state. The analytical development in 
the present case follows the same procedure as in the previous sec- 
tion. However, it is more laborious, due to the presence of an ad- 
ditional parameter which is the phase (f>{r) of the periodic equi- 
librium state. We shall refer to reference [4] for complete de- 


tails. We give the result obtained for the lift force 


° L (t) = H'/rC? 


i iT( Tr +a Ta ] W 


f dra(r)C[ r ' ln, |tt( r),a(i 

Jo 


where C e ^ mi denotes the equilibrium component of the lift force and 
CV ant denotes its transient part. In the case of vortex- induced 
vibrations, the value of C£ JU,i can be deduced from experimental data 
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on a stationary cylinder: 

C e L qml = yl 00 :a(T)]sin{A: s jo:(r)](/ - r) ^[a(r), a(r), r]} (45) 

with yloo and fcs known from experimental data. 

The contribution into Ci(t) contains terms associated with C‘ L qml (the 
first two terms of the right-hand-side member of equation (44)). 
Their contribution may exceed the value of C e ^ ml (a(t), <j>(t);t, t = t ) . 
as shown by our own numerical simulation, based on the following 
expression for <f>(r) : 

<p(r) = Apa(r ) + Bpa(r) (46) 

Such relation for <j>{~) is suggested by physical reasoning [4] . 

It appears clearly that relation (43) , giving the lift force in the 
first physical situation with time-invariant equilibrium state, is 
not adequate for the periodic equlibrium state. Modeling of sys- 
tems in flight dyanmics involving a new regime (for instance, stall 
phenomenon) requires an appropriate change of the analytical ex- 
pression of the lift force. 

5. Predictions for a periodically forced cylinder near the Hopf 
bifurcation: relationship with the amplitude equation approach. 
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It appears from relation (44) the the lift force depends generally 
on the motion history of angle of attack, which is itself depen- 
dent on the value of the lift force (through the governing equa- 
tion of motion of the elastically mounted cylinder). So, gener- 
ally the equations governing the angle of attack and the lift force 
are coupled. In the special case of a periodically forced cylin- 
der, one controls externally the motion and therefore the angle - 
of -attack. Hence, the only independent variable of the system is 
constituted by the lift force, i.e. the velocity field. One can 
simplify furthermore the analysis by restricting it just near the 
Hopf bifurcation. Such a restriction corresponds to the hypoth- 
esis of the amplitude equation approach. 

The equations describing the present physical situation correspond 
to equations (26) slightly modified: the terms Cj.(l — H{t+ )) are 

to be replaced by appropriate terms associated with external forc- 
ing (cf. [4]). Due to the restriction of the analysis to the neigh- 
borhood of the critical situation, only the equation correspond- 
ing to the projection onto the critical eigenvalue 70 (associated 
with the eigenvalue u/ 0 ) is interesting to study. We make use of 
the theoretical results of Elphick et al. [6] to carry out the sim- 
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plification of this equation to get its ’‘normal form’’ (cf. ref- 
erence [4] for complete details). The final equation obtained is: 
dA , , 

— = {a + two) A - b\A\*A -r cf 0 e' ut (47) 

where a, 6 are complex constants, c is a real constant, f 0 is the 
forcing amplitude and u: is the external frequency. 

It can be shown that equations (47) are derived from the follow- 
ing modified Van der Pol equation: 

A - (a - jA 2 )A + u> 2 A - rjA 3 = / 0 sin (wf) (48) 

Such an equivalence was already noticed by Provansal [7]. However, 
the author hadn’t realized its relationship with oscillator mod- 
els of vortex- induced oscillations. Indeed, one can show that the 
lift force generated by vortex- shedding behaves as the amplitude 
of the velocity field. 

We have therefore justified the choice of Van der Pol equation as 
the model of the lift force generated by vortex- shedding from fluid 
dynamics considerations (Havier-Stokes equations). Furthermore, 
our analysis supports the analytical form of the fluid oscillator 
proposed by Skop and Griffin [8] . We have also analysed the am- 
plitude equation (47) to show that there is jump and hysteresis phe- 
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nomena inside the lock-in regime (where the body oscillates at the 
forcing frequency) . 

6 . Conclusion 

In the present study, we have taken some effort to improve our math 
ematical modeling approach in order to incorporate a Hopf bifur- 
cation. We have developed the indicial response of the velocity 
field, based on the incompressible Navier-Stokes equations for the 
regimes of time -invariant and periodic equilibrium state. 

The theory has been applied to the particular case of flow past a 
cylinder in periodically forced motion. It has been shown that the 
indicial response approach is equivalent to the currently popular 
amplitude equation approach, near the Hopf bifurcation. By rec- 
ognizing that the lift force behaves like the amplitude of the ve- 
locity field, we have been able to justify the choice of a Van-der- 
Pol equation which was used empirically to model the lift force gen 
erated by vortex-shedding. We have been also able to capture jump 
and hysteresis phenomena that experiments indicate occur within the 
lock-in regime (where the velocity field oscillates at the forc- 
ing frequency) . 

It is hoped that our modeling approach will have some beneficial 
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implications for mathematical modeling of other configurations of 
self -induced vibrations that bear common features. 
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LEGEND . 


Fig.l.- Maneuver referred to space-fixed {X,Y) and moving (x, y) 
ordinates, the latter attached to the airplane. 

Fig. 2.- Formation of indicial response. 
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